Final Mass and Maximum Spin of Merged Black Holes and the Golden Black Hole 
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We present results on the mass and spin of the final black hole from mergers of equal mass, 
spinning black holes. The study extends over a broad range of initial orbital configurations, from 
direct plunges to quasi-circular inspirals to more energetic orbits (generalizations of Newtonian 
elliptical orbits). It provides a comprehensive search of those configurations that maximize the 
final spin of the remnant black hole. We estimate that the final spin can reach a maximum spin 
a/Mh ~ 0.99 ± 0.01 for extremal black hole mergers. In addition, we find that, as one increases 
the orbital angular momentum from small values, the mergers produce black holes with mass and 
spin parameters {Mh/M, a/Mh} spiraling around the values {Mh/M, a/Mh} of a golden black hole. 
Specifically, (M^ — Mh)/M ex. e^^"^ cos(f> and (a — a)/Mh oc e^'^'^sin^, with cp a monotonically 
growing function of the initial orbital angular momentum. We find that the values of the parameters 
for the golden black hole are those of the final black hole obtained from the merger of a binary with 
the corresponding spinning black holes in a quasi-circular inspiral. 

PACS numbers: 



Introduction: An important question that can only be 
answered with the tools provided by numerical relativ- 
ity is: What is the final mass and maximum spin of the 
final black hole (BH) from generic binary mergers? In 
Ref. fH , referred to here as Paper I, we took the first step 
towards answering this question. We studied the behav- 
ior of a 1— parameter series of equal mass, non-spinning 
BHs in merger or fly-by. The adjustable parameter was 
the magnitude of initial linear momentum P for each of 
the BHs. We kept fixed the initial binary coordinate sep- 
aration d = 10 M X and the angle 9 = 26.565° between P 
and d. That is, the binaries had initial orbital momen- 
tum L — dP sin^z — AA7 M P z and impact parameter 
b = d am9 = AAl M. Here M = 2 m is the total mass of 
the binary, where m is the mass of each BH. The main 
result in Paper I was the specific dependence or trans- 
fer function we found that connects the final merged BH 
spin a/Mh with the initial L/M^^^ (see Fig. 1 in Paper 
I, which is repeated in Fig. [T] here in the curve labeled as 
S/]\/P — 0.0, i.e. vanishing total initial spin). For this 
case, a maximum final BH spin a/Mh ~ 0.82 was found 
at L/Ml,„^ « 1. 

The result of Paper I has particular astrophysical in- 
terest as an upper limit on the spin from initially non- 
spinning BH mergers, though BHs can be spun up by 
other processes as well, specifically by accretion. No 
previous BH merger simulation has produced final spins 
close to maximal. We do know from specific examples 
in work by Campanelli et al. \2\ that "hangups" occur, 
in which the orbital radius decreases only slowly or not 
at all, and which enable the radiation of substantial an- 
gular momentum before the merger. In Paper I, we also 
found a splash-skip behavior involving what Pretorius 
and Khurana H called whirl orbits. In these cases the 



interacting BHs can approach, then recede to a large dis- 
tance, then approach again, producing repeated bursts of 
gravitational radiation. This is clearly an extreme case 
of "hangup." 

Here we extend the work in Paper I to encounters 
where at least one of the BHs carries a spin, with the 
spins parallel or anti-parallel to the orbital angular mo- 
mentum. The main motivation is to add extra initial an- 
gular momentum in order to increase the final BH spin. 
In addition, we vary the angle 0, but keep the initial sep- 
aration as before, d = 10 M. We find similar transfer 
functions to that in Paper I, with the maximum final BH 
spin depending on the spins of the merging BHs and the 
scattering angle 9. Our simulations show that the final 
BH can reach a maximum spin a/Mh ~ 0.99 ±0.01 when 
extrapolated to the merger of maximally spinning BHs. 

In addition, we found an interesting behavior regarding 
the parameters {Mh/M, a/Mh} of the final BH as we 
increased the orbital angular momentum. The values for 
{Mh/M, a/Mh} spiral around the parameters of what we 
operationally call a golden BH. We have been able to 
identify the golden BH to be the hole obtained from the 
merger of a binary in a quasi-circular inspiral with the 
corresponding spinning BHs. 

The computational setup for the simulations presented 
here is very similar to the one in Paper I, the main dif- 
ference being the spins of the colliding BHs. All sim- 
ulations had 10 levels of refinement, sixth order spatial 
differencing, with an outer boundary of ^ 320 M. For 
the runs used to compute the 9 > 26.565 maximum final 
spin 200 — 300 M runtimes) a resolution of M/51.6 
was used on the finest level, with each subsequent level 
decreased by a factor of 2. For the 9 < 26.565° runs 
(~ 200 - 300 M runtimes), a resolution of M/103 was 
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used. For the longer runs to compute the golden BH 
500 — 3000 A/ runtimes), a resolution of A//90 was 
used on the finest grid. 

The present study required f , 000 binary black holes 
(BBH) simulations. Every data point displayed in each 
figure involves at least one BBH simulation and in Figs. [2] 
and O up to ten simulations. It was thus prohibitively 
expensive to do the runs needed to estimate errors via 
Richardson extrapolation for each data point. With this 
in mind, we chose to do all the simulations at the reso- 
lutions mentioned in the previous paragraph, which by 
experience we know our code yields sufficiently accurate 
results. In order to get error estimates, we selected a 
subset of representative cases, in some instances the most 
challenging cases (e.g. high BH spins), to perform conver- 
gence tests and apply Richardson extrapolation. Those 
are the cases with error bars in Figs. [1] |3] and [SJ The 
convergence tests demonstrated consistency with the ex- 
pected fourth order accuracy of our code. In addition, the 
Richardson extrapolation of results yielded error bars for 
the spin and mass of the final black hole ~ 0.1%. For 
situations involving BH spins a/Mh > 0.95, we identified 
additional sources of error due to inaccuracies in find- 
ing the exact location of the apparent horizon (AH) that 
yield errors in the range between 0.5 and 1%. 

Maximum Final Black Hole Spin: The first compo- 
nent of this work is aimed at addressing the following 
three questions regarding the maximum spin of the fi- 
nal BH: Does the maximum occur near the same initial 
angular momentum as found in Paper I? What is the de- 
pendence of the maximum residual spin on the value of 
the initial spins of the merging BHs? How close is the 
maximum final spin to the Kerr limit? The answer to 
the first question is "yes" . The maximum for the resid- 
ual angular momentum occurs near J/M^^^ w 1 for 
all cases considered. This conclusion is apparent from 
Fig. [TJ where we show the dependence of the spin of 
the final BH a/Mh as a function of the initial total (or- 
bital plus spin) angular momentum J/M^^^ at a con- 
stant 6 — 26.565°. From top to bottom are the cases 
S/]VP = {0.4,0.3,0.2,0.1,0.0,-0.1,-0.2,-0.3,-0.4}, 
respectively, where S/A'P = (mai + ma2)/(2m)^ = 
(fli/m + a2/m)/A is the total initial spin of the binary 
with fli/m and a2/m the individual BH dimensionless 
spin parameters. For reference, in the S/A-P — 0.4 case, 
we show error bar estimates obtained from Richardson 
extrapolation. 

An interesting feature in the results displayed in Fig.[T] 
is the crossing of lines for negative S/M"^. The cross- 
ing occurs because the mergers with negative S/M"^ do 
not suffer the "hangup" of the positive counterparts . 
Therefore, as S/M"^ becomes more negative, one is able 
to increase the total angular momentum J/M'^^^ to 
higher values before reaching the maximum spin ajMy^ 
which roughly corresponds to the transition between di- 
rect plunges and inspiral-like coalescences. 

We have looked for but find very little torque-up, sug- 
gesting very little spin-orbit coupling. At the same time. 




FIG. 1: Final spin a/Mh vs the initial total (orbital plus 
spin) angular momentum J/Madm^ at a constant 9 = 
26.565°. From top to bottom are the cases S/M"^ = 
{0.4, 0.3, 0.2, 0.1, 0.0, -0.1, -0.2, -0.3, -0.4}, respectively. As 
a reference, error bar estimates are given for the S/M^ = 0.4 
case. 



for the configurations we consider, the final spin as a 
function of J depends only on the total initial spin, 
not how it is distributed between the interacting BHs. 
For instance, the cases 01,2/^ — {0.4,0.4} {S/M"^ = 
0.1 + 0.1), ai,2/m {0.0, 6.8} {S/M^ = 0.0 + 0.2) and 
01^2/ m = {0.6, 0.2} (S'/A/2 = 0.15+0.05), which have the 
same initial total spin S/AI^ — 0.2 but distributed dif- 
ferently between the two BHs, produce almost identical 
final merged BH angular momentum vs. the initial total 
angular momentum J/M^^^. We have carried out sim- 
ilar experiments for total S/M'^ = 0.1 (with a\_2/m = 
{0.8,-0.4} and ai.2/m = {0.2,0.2}) and S/M^ = 
(with ai,2/m = {6.0,0.0}, ai,2/m = {0.4,-0.4} and 
01,2/"^ — {0.8,-0.8}). These results suggest that the 
initial spin is simply "bundled into" the total J in the 
evolution. A fraction of J is then deposited into the final 
BH, so that the final BH spin depends essentially only 
on the total initial angular momentum. 

This result is confounding because a close examina- 
tion of the orbits shows significant differences between 
two cases with the same total but different distribution 
of initial spin. For instance, the case ai^2/'m — {0.4, 0.4} 
with symmetrical initial data yields a significantly differ- 
ent evolution from the case ai,2/m = {0.0,0.8}. In fact, 
the a\;2/m = {0.0,0.8} case produces a kick of ^ 175 
km/sec, and the symmetric case 01^2/"^ = {0.4, 0.4} pro- 
duces no kick 0-0|- 

We argued in Paper I that the existence of the max- 
imum in the final spin parameter a/Mh depends on the 
participation of an "intermediate excited state" which we 
characterized as essentially a highly distorted BH. Exis- 
tence of such an intermediate state seems to be generic, 
and it emits the largest part of the radiated energy 
and angular momentum. We conjecture that the non- 
linearity of such a state will produce such strong radia- 
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tion of angular momentum that it will enforce an upper 
limit on the residual BH angular momentum. We extend 
that argument here to claim that, while in this state, sub- 
stantial information is lost about the initial configuration 
of the system. 

The answers to the second question ( "how does the fi- 
nal BH depend on the initial total spin?" ) can be found 
in Fig. [5] where we show the maximum spin a/Mh de- 
posited in the final BH vs the initial total spin angu- 
lar momentum S/M"^. From top to bottom are 9 — 
{26.565°, 32°, 45°, 55°, 70°, 80°, 90°}, respectively. Each 
9 = constant line in Fig. [5] is obtained from the set 
of maximum final spins like those found in Fig. [1] for 
9 = 26.565°. Notice in Fig. [5]that the dependence of the 
final spin on the spins of the merging BHs is almost linear 
for 9 approaching 90° but non-trivial for small values of 
9. The small 9 values correspond to small impact param- 
eters. For these cases, mergers with negative S/A'P are 
direct plunges; very little angular momentum is radiated, 
thus the turn around observed in Fig. [2j 

A possible answer to the third question ("what is the 
maximum final spin in these mergers?") can be found 
by extrapolating the 9 = constant lines in Fig. [2] to max- 
imally spinning incident BHs (i.e. to S/M"^ = 0.5). Such 
extrapolation yields a final BH spin a/Mh ~ 0.98. This 
maximum final spin is slightly lower than the predic- 
tion for quasi-circular inspirals by Kesden 8] of a/Mh ~ 
0.9988 but higher than the estimate by RezzoUa W . of 
a/Mh ~ 0.959. Recently, work by Sperhake et al. 
has also produced final spins a/Mh ~ 0.95 for non- 
spinning, highly boosted BHs. We believe that these 
highly boosted BHs address a different regime from the 
one considered here. The mergers in Sperhake et al. [lol | 
involve ultra-relativistic encounters with initial total an- 
gular momentum as high as J/M^^^^ « 3 and small 
enough impact parameter to "force" a merger; that is, 
the binaries are initially "over-saturated" with angular 
momentum. In our case, we start with J /M'^^^ k, 1 
configurations and arrange the binary to minimize the 
radiation losses of angular momentum. 

Another important aspect to keep in mind regarding 
the search for the maximum final spin is that our study 
involves a 2-parameter family of simulations in S/M"^ and 
9. Extrapolation along 9 = constant as in Fig.[2]may not 
be the right path in parameter space to find the maxi- 
mum final spin. To further investigate this point, we plot 
in Fig. [3] the final spin as a function of 9 but now for lines 
of S/M"^ = constant. From bottom to top are the cases 
S/M'^ = 0.0,0.1,0.2,0.3,0.4. The crosses along the top 
axis of the figure are the values obtained from the 9 = 
constant extrapolation to extremal BHs in Fig. [2] It is 
now clear why the maximum final spin values obtained 
from extrapolating the data in Fig. [5] are similar to each 
other (i.e. all lines approach each other for S/M"^ = 0.5). 
As Fig. [3] shows, the final spin for 9 > 30° along S/M"^ = 
constant lines is roughly constant thus yielding similar 
extrapolation values along 9 — constant (see crosses in 
Fig. El). 
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FIG. 2: Maximum final spin vs the initial spin angular mo- 
mentum S/M^. From top to bottom are the cases 6 — 
{26.565° , 32° , 45° , 55° , 70° , 80° , 90° }, respectively. 




FIG. 3: Maximum final spin vs 9 for S/M^ = 
0.0, 0.1, 0.2, 0.3, 0.4 (bottom to top) 



Notice that in Fig. [3] we have extended the range of 
simulations to include 9 < 26.565°. A completely differ- 
ent 9 dependence is found for these smaller angles, whose 
details are better appreciated in the inset of Fig. [31 where 
we have also included the error bars estimates for the 
S/M^ = 0.4 case. As 9 decreases along each S/M^ = 
constant line, the final spin increases and reaches a max- 
imum. Notice in particular from the inset Fig. |31 that 
the top line, S/M'^ = 0.4 case, is not as smooth as the 
others due to an artificial loss of angular momentum like 
that seen by Marronetti et al. [ll|. We however are able 
to obtain a maximum spin of a/Mh = 0.98 ± 0.01 around 
9 — 16°. Notice then from Fig. [3l that a better estimate 
of the maximum (and hopefully global) final spin can be 
found by an extrapolation to extremal BHs using the val- 
ues of the maximum on each the S/M"^ = constant lines 
in Fig. [3l The result of this extrapolation yields a max- 
imum final spin a/Mh ~ 0.99 ± 0.01 at 6* ~ 15.2° and is 



4 




0.892 



1.3 1.35 



FIG. 4: Final mass Mh/M (top panel) and spin a/Mh (bot- 
tom panel) vs JjM^^^ for Q = 60° (solid), 80° (dash) and 
90° (dotted) with initial total spin of SjU"^ = 0.4. The 80° 
and 90° results almost exactly overlap. 
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FIG. 5: Final mass Mh/M (dashed, top) and spin a/Mh (dot- 
ted, bottom) vs J/M^^^ with error bars for 9 — 90° and 
S/M^ = 0.4. 



depicted as a black dot in Fig. |3l Also very important 
to point out is that, although this study involved simula- 
tions in the highly non-linear regime of general relativity, 
the data used to extrapolate to maximum final spin ex- 
hibit only a slight departure from linearity (see maximum 
values in the inset of Fig.[3]and upper right corner data in 
Fig. [2]). In summary, our study consists of a 2-parameter 
{S'/M^, 6} family of simulations, with global maximum 
a/Mh ~ 0.99 ± 0.01 of the final spin found by slicing the 
data along S/K'P — constant lines. 

Golden Black Hole: In addition to studying the max- 
imum possible final spin, we investigated the final state 
{Mh/M, a/Mh] as one increases the initial angular mo- 
mentum J/Ml^^. The results for the case S'/A/^ 0.4 
with (ai,2/m = 0.8) and angles 9 = {60°, 80°, 90°} can 
be found in Fig. S) In order to get a sense of the accu- 
racy of our results, we repeat in Fig. [5] the case 9 — 90° 
and S/M'^ = 0.4 including error bar estimates for both 



the final mass and spin. It is interesting that as J /M"^^^ 
is increased, the pair {Mh/M, a/Mh] exhibits decaying 
oscillations around some values {Mh/M, a/Mh]- After 
a close examination of the orbits that lead to these re- 
sults, we concluded that these oscillations are a direct 
consequence of the binary having to undergo multiple 
zoom-whirl episodes to shed the "excess" of angular mo- 
mentum that inhibits the merger [l^. Note also that 
for a given angle 9, the amplitude of the oscillations are 
closely similar between a/Mh and Mh/M . The two quan- 
tities are 90° out of phase. Furthermore, the frequency 
of the oscillations is not constant but increases for larger 

In order to get a better insight on the oscillation and, 
in particular, to investigate whether the damping of the 
oscillations continues for larger J/M^^^, we focused our 
attention on the case = 0. Vanishing initial spins 

not only yield shorter merger times but also require lower 
resolutions, thus reducing the computational cost of the 
large number of simulations needed. Figure [6] shows 
Mh/M (top) and a/Mh (bottom) vs J/Ml^^ for 9 = 90° 
and non-spinning initial BHs. Crosses are the data from 
the simulations and solid lines are fits to 
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with (p a monotonically growing third order polynomial 
in J/M^j^^. The fittings were done as follows: Around 
J*IMldm ~ 0-985 (vertical line in Fig. E]) the oscilla- 
tion is almost completely damped out. We have then 
divided the data into values below and above J*/Af^^„j. 
For those below, we fit a decaying exponential and, sim- 
ilarly, a growing exponential for the data above values 
J*/M^^^. For the fitting, we have ignored the four points 
in the neighborhood of J*/A/^^^ since in that region the 
oscillations are almost completely damped. 

An interesting finding is that, modulo a phase shift, (p 
is basically the same fitting function (i.e. similar constant 
coefficients in the third order polynomial in J/M"^^^ 
for the decaying and growing sectors. We also obtain 
that both the decaying and growing fittings yield the 
same values of {Mh/M, a/Mh] = {0.951,0.685}. The 
only changes are in the constants {A±, B±, C±,D±]. Fi- 
nally, the exponential envelopes (dashed lines in Fig. |6]) 
cross approximately at J*/M'^^^. For this non-spinning 
situation, we repeated similar analysis for 9 = 70° and 
found the same values for {Mh/M, a/Mh] within our nu- 
merical errors. In addition, we considered total initial 
spins cases S/M"^ = 0.2 {9 = 90°) and S/M^ = 0.4 
{9 = 60°, 80°, 90°). The S/M'^ = 0.2 case yielded 
{Mh/M^/Mh] = {0.936,0.800}, and the S/M^ = 0.4 
values{M?,/M,a/M,,} = {0.906 ± 0.001,0.900 ± 0.001}, 
with the last case selected for^ror estimates. We name 
the BHs with parameters {Mh/M,a/A'Ih] the golden 
BHs. 
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FIG. 6: Mh/M (top) and a/Mu (bottom) vs J/Ml^^ for 
e = 90° with initial total spin of S/M^ = 0.0. 




Mj/M 



Figure [7] combines the results from Fig. [6] The vertical 
line denotes the location of the golden BH. The result is 
a logarithmic spiral which approaches the golden BH as 
'V^^ldm nears ■h/M'^^^^ « 0.985 from above or below. A 
similar spiral is found in the a/Mh vs M^/M plane when 
the angular momentum J/M^^^^ is held constant and one 
plots them instead as a function of the angle 6. The cen- 
ter of the spiral is again at the golden BH. A closer look 
at our simulations reveals the identity of the golden BH 
as the final BH obtained from the merger of a binary of 
the corresponding spinning BHs in a quasi-circular inspi- 
ral. We are currently investigating the reasons behind 
this finding. 

Conclusions: From the first studies of binary BH merg- 
ers jlSrtlS], it became apparent that merger waveforms 
have a simpler structure than was expected. Here we 
have presented two more examples of reduction of com- 
plexity. One is the invariance of the final BH spin vs the 
initial total angular momentum with respect to how the 
individual spins are distributed betweenjihe interacting 
BHs. The other invariance is the final {Mh/M, a/Mh} of 
the golden BH with respect to the angle 9, for a given 
total initial spin S/M"^ (see Fig. . These are explicit 
reductions of degrees of freedom from the initial data to 
the final BH. From the results of our 2-parameter fam- 
ily of simulations, we estimated that the maximum final 
spin is found for initial extremal BHs with spins aligned 
with the orbital angular momentum and 6 w 15.2°. Our 
results imply a maximum final spin a/Mh ~ 0.99 ± 0.01. 
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